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The study of the Extragalactic Background Light (EBL) is crucial to understand many astro- 
physical problems (as the formation of first stars, the evolution of galaxies and the role of dust 
emission). At present, one of the most powerful ways to put constraints on EBL is represented by 
the study of the photon-photon absorption on gamma-ray spectra of TeV blazars. Adopting this 
method, we found that, if the only contribution to the optical and Near Infrared (NIR) background 
is given by galaxies, the spectrum of the blazar H 1426+428 cannot be fitted. To reproduce the ob- 
servational data of H 1426+428 a Near Infrared excess with respect to galaxy counts is required, 
with amplitude consistent with both the Matsumoto et al. (2000) data with Kelsall's model of 
zodiacal light (ZL) subtraction and the DIRBE data with Wright's model of ZL subtraction. The 
derived constraints on the optical EBL are weaker, because the experimental errors on blazar data 
are still bigger than the differences among various optical EBL models. In the mid-infrared the 
SPITZER measurement of v/ v =2.7 nW raT 2 sr~' at 24 /jm, provides the best fit of the blazar 
spectrum. 



Baryons in Dark Matter Halos 
5-9 October 2004 
Novigrad, Croatia 



* Speaker. 

' We thank A. Celotti, L. Costamante, S. Inoue and T. Matsumoto for useful discussions. 



Published by SISSA 



http : / /pos .sissa.it/ 



Extragalactic Background Light and photon-photon absorption 



Michela Mapelli 



1. Introduction 

The study of the Extragalactic Background Light (hereafter EBL) might provide unique in- 
formation to understand many crucial astrophysical questions, including, among others, the first 
cosmic star formation (Salvaterra & Ferrara 2003), the evolution of galaxies and the role of dust 
emission (Totani & Takeuchi 2002). In this proceeding we want to show as the study of the photon- 
photon absorption on blazar spectra can provide a powerful method to improve our knowledge of 
the EBL. In Sec. 2 we summarize the present measurements and problems about EBL. In Sec. 3 we 
shortly describe the photon-photon absorption and our procedure. In Sec. 4 we present our findings 
and in Sec. 5 we report some conclusions. 

2. The Extragalactic Background Light 
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Figure 1: Extragalactic Background Light data and corresponding models. Left panel: MK1 (dashed blue 
line), MK2 (solid blue line). The Matsumoto et al. (2000) data are represented with filled triangles.Right 
panel: DW1 (dashed red line), DW2 (solid red line). The Wright (2001) data are represented with open 
circles. In both the two panels are shown: CI (dotted green line), data from Bernstein et al. (filled squares), 
Madau & Pozzetti (filled circles), Metcalfe et al. (cross) and Papovich et al (open triangle). 

At present, our knowledge of the EBL still presents a lot of uncertainties, both in the optical 
and in the infrared region. Galaxy counts (Madau & Pozzetti 2000) could entirely represent the 
optical EBL. But Dube et al. (1979) derived an upper limit of the background flux at X ~ 5000A 
considerably higher than galaxy counts and, more recently, Bernstein et al. (2002) measured a 
significant excess at X=3000, 5500, and 8000 A (Fig. 1). On the other hand, this excess can not yet 
be confirmed, as pointed out by Mattila (2003). An analogue uncertainty about the presence of an 
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excess affects the Near-Infrared region (hereafter NIR). The principal problem of NIR background 
measurements is represented by the subtraction of Zodiacal Light (hereafter ZL). The Zodiacal 
Light is the sunlight scattered by interplanetary dust. There are two different models of ZL light 
subtraction: the Kelsall model (Kelsall et al. 1998), based on the temporal variability of the ZL due 
to the different amounts of dust encountered by the Earth during its orbit, and the Wright model 
(Wright 2001 and references here), the maximal one, based on the fact that the Kelsall model seems 
to be unable to subtract all the ZL at some wavelengths (for example at 25 /urn). In the left panel 
of Fig. 1 we indicate the Matsumoto data (Matsumoto et al. 2000) of the NIR region, obtained 
subtracting the ZL contribution through the Kelsall model; while in the right panel we show the 
COBE/DIRBE data (Wright 2001), where the ZL has been subtracted following the Wright model. 
Both Matsumoto et al. and Wright data show an excess in the NIR with respect to galaxy counts, 
even if of different amount. 

The problem of the Mid and Far-Infrared region is the small number of data. In fact the only recent 
data are provided by ISOCAM (Metcalfe et al. 2003), measuring a flux of 2.7 ± 0.62 nW m~ 2 sr _1 
at 15/an, and by SPITZER satellite (Papovich et al. 2004), which gives a background total flux of 
2-7^Qg nW m 2 sr~ 1 at 24 /mi. At present we are waiting for new SPITZER measurements. 



3. The photon-photon absorption on blazar spectra 

Because of the present uncertainties about EBL direct measurements, we need an alternative 
approach to improve our knowledge of the Extragalactic Background Light. This alternative ap- 
proach can be represented by the study of the photon-photon absorption on blazar spectra, a pair 
production process involving two photons whose energies, E and £, are related by the (approxi- 
mate) formula £ ~ 2(m e c ) /E (where m e is the electron rest mass). This means that, if E is the 
energy of a photon in the TeV range of a blazar spectrum, the interacting photon should belong to 
the Infrared Background (X ~ 2(£'/lTeV) i um), where X is the wavelength of the photon of energy 
£). This strong relation between the energies of the interacting photons exists because the photon- 
photon absorption cross section for a given E is a well peaked function of X (Fig. 2). The full 
expression of the photon-photon absorption can be written (Madau & Phinney 1996) as 

r Zem dl r [ ( l — x) f°° 

i(E)= I dz— dx- -/ den(e)o(e,E,x), (3.1) 

Jo dzJ-i 2 Je lh 

where dl/dz is the proper line element 1 jj| = [(1 +z) £ (z)]~ , c is the speed of light and 

E(z) = [£l M (l+z) 3 +n A +(l-n A -n M )(l+z) 2 ] 1/2 - Ineq. |3l[ x = cos 6, 6 being the angle 
between the directions of the two interacting photons. E = Eq (1 +z) is the observed energy of the 
blazar photon and £ = £o(l +z) is the observed energy of the background photon; z em is the redshift 
of the considered blazar; finally, n(z) is the specific number density of background photons. The 

2 m 2 c 4 

energy threshold for the interaction, £,/,, is defined by £ ; /, = E ^{_ x ^ . where m e is the electron mass. 
The photon-photon absorption cross section is given by 



a(£,£,x) = ^(l-p 2 



2p(p / -2) + (3-p 4 )ln 



(3.2) 



,1-P, 

'We adopt the following cosmological parameters: Hubble constant Hq=11 km s Mpc~', £2^=0.27, £"2^=0.73, 
which are in agreement with the recent WMAP determination (Spergel et al. 2003). 



/3 



Extragalactic Background Light and photon-photon absorption 



Michela Mapelli 




Figure 2: Photon-photon absorption cross section as a function of the background photon wavelength, A, 
integrated over the interaction angle and over the redshift (with Zem = 0.129, the redshift of H1426+428), 
for three values of the observed energy of the blazar photon: 100 GeV (dashed line), 1 TeV (dotted) and 10 
TeV (solid)). 



where Oj is the Thompson cross section and p 
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We integrated numerically the 



expression [3JJ for different EBL models, using a three-dimensional integral based on the method 
of Gaussian quadratures. 

The absorbed spectrum of the blazar ({j§) abs ) can be derived by the unabsorbed one ( b ), 
given the relation: 

'dN\ 



dE J 



-i(E) 



abs 



dEJ 



(3.3) 



unabs 



Assuming that the unabsorbed spectrum has a power law shape 2 ((dN /dE) una b s °< E~ a ) and cal- 
culating x(E) for various EBL models, we can derive the absorbed spectrum and compare it with 
observational data. Through this procedure we are able to distinguish among various EBL models. 
We apply this method to 5 different EBL models, reported in Table 1. 



4. Results 

The first result we obtained adopting our method is that an EBL excess (in the optical or NIR) 
is required to fit the TeV spectrum of the blazar H 1426+428 (redshift z = 0.129). In fact adopting 

2 The assumption that the unabsorbed spectrum follows a power law is valid especially in the case of the blazar we 
are going to analyze (H 1426+428). For other blazars the unabsorbed spectrum is better fitted by a power law with an 
exponential cut-off (Konopelko et al. 2003). 
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Model Optical Background 



NIRB 



MIRB-FIRB 



CI Madau & Pozzetti (2000) 

MK1 Madau & Pozzetti (2000) Matsumoto et al. (2000) (K) a 

MK2 Bernstein et al. (2002) Matsumoto et al. (2000) (K) 

DW1 Madau & Pozzetti (2000) Wright (2001) (W) a 

DW2 Bernstein et al. (2002) Wright (2001) (W) 



Totani & Takeuchi (2002) + SPITZER h 
Totani & Takeuchi (2002) + SPITZER 
Totani & Takeuchi (2002) + SPITZER 
Totani & Takeuchi (2002) + SPITZER 
Totani & Takeuchi (2002) + SPITZER 



Table 1: Summary of the considered EBL models. 



"(K) and (W) indicate the ZL subtraction obtained using Kelsall's model and Wright's model, respectively. 
^Totani & Takeuchi model rescaled to the SPITZER data (Papovich et al. 2004). 



the model CI, i.e. assuming that both the optical and the NIR background are represented by galaxy 
counts, we cannot fit the spectrum of H 1426+428, as you can see from Fig. 3 and from the % 2 test 
in Table 2. 
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Figure 3: Three different fits to the specttum of H 1426+428. The photon-photon absorption has been 
calculated assuming the background model CI (see Table 1). The shaded area indicates the region between 
the absorbed spectra obtained with the upper limit (+lo) of galaxy counts and that obtained with the lower 
limit (-lo) of galaxy counts. The solid line represents the case CI (see Table 1). The observational data 
reported here are from CAT 1998-2000 (filled squares), Whipple 2001 (filled triangles), HEGRA 1999-2000 
(filled circles), HEGRA 2002 (open circles). All error bars are at 1 a (Aharonian et al. 2003). 

Then we need of an EBL excess. If we assume that this excess is in the NIR region, we can 
fit the spectrum of H 1426+428 (Fig. 4), both if we calculate %(E) adopting the Kelsall model of 
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CI MK1 


MK2 


DW1 


DW2 


X 2 48.00 18.69 


19.71 


6.98 


7.41 



Table 2: % 2 values of the considered models. 
We calculated the % 2 for each one of the considered models. We based our statistical analysis on 12 observational data 
of the spectrum of H1426+428 (CAT 1998-2000, Whipple 2001, HEGRA 1999-2000). 

ZL subtraction for the NIR (MK1 and MK2 models) and if we calculate x(E) adopting the Wright 
model (DW1 and DW2). The DW1 model gives the best result, considering the % 2 (Table 2); but 
we cannot exclude also the other 3 cases (DW2, MK1 and MK2). 

Unfortunately, our results are not conclusive concerning the existence of an optical excess. The 
reason is, as shown in Fig. 4, that the difference between the results obtained using only galaxy 
counts for the optical background and those obtained adopting the Bernstein et al. measurements 
is much smaller than the error bars in the observational data of the blazar spectrum. 




0.1 1 0.1 1 10 

E 7 (TeV) 

Figure 4: Fit to the H 1426+428 spectrum. The EBL models assumed are: Left panel: MK1 (dashed blue 
line) and MK2 {solid); Right panel: DW1 {dashed red line) and DW2 {solid). Shaded areas refer to ±lo 
errors, the observational data are the same as in Fig. 3. 

Using our method we can also derive the most likely spectral index of the blazar spectrum. In 
particular, if we use the DW1 model, we obtain for H 1426+428 a spectral index a = 1.8, in good 
agreement with previous results (Aharonian et al. 2003). 

Our method can provide also some constraints on the SPITZER data, at least under the assump- 
tion that the Totani & Takeuchi (2002) model of Mid and Far Infrared is correct. Let's consider 
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the SPITZER measurement at 24 ^m (total background flux F = 2.7^ ' 7 nW m _2 sr _1 ). We obtain 
a good fit of the blazar spectrum only adopting the SPITZER best value; while using the upper 
or the lower limit of this measurement our results are not in agreement with some observational 
data (Fig. 5). 
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Figure 5: Fit of the H 1426+428 spectrum. We have assumed an EBL model given by: (a) optical back- 
ground (0.3-1.2 pm): Madau & Pozzetti 2000; (b) near infrared background (1.2-4. pm): DIRBE (Wright 
2001) with the ZL model of Wright & Reese (2000); (c) middle and far infrared background : model of 
Totani & Takeuchi (2002). The model of Totani & Takeuchi for the Mid and the Far Infrared has been 
rescaled in the 8-30 pm range assuming different values of the EBL at 24 pm, in particular assuming: 2.7 
(red line), 3.8 (blue line) and 2 (green line) nW m~ 2 sr -1 . The observational data reported here are from 
CAT 1998-2000 (filled squares), Whipple 2001 (filled triangles), HEGRA 1999-2000 (filled circles). 



5. Summary 

In conclusion, the most important finding of our work is that a Near Infrared excess with 
respect to galaxy counts must be present, if we want to explain the photon-photon absorption 
deduced from blazar spectra. We cannot definitely distinguish between the NIR excess derived 
adopting the Kelsall model of ZL subtraction and that obtained from the Wright model; but the % 2 
analysis favors the latter. This result is crucial for models describing the connection between EBL 
and Population III stars (Salvaterra & Ferrara 2003). In fact these models, which explain the EBL 
in the Near Infrared region as the radiation emitted by high z metal free stars and redshifted in the 
infrared region, are in agreement with the existence of a NIR excess with respect to galaxy counts. 

Our work presents also a lot of challenges. First of all we need new measurements in the 
0.1-1 TeV range of the spectrum of H 1426+428, because this region is crucial to distinguish 
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among various optical background models. These new measurements could be provided by new 
Cherenkov telescopes, which are already at work (MAGIC, VERITAS, HESS) or will be 
operating in a few years (GLAST). We need also new spectra of other high redshift blazars, 
because now our analysis is based on only one blazar (H 1426+428). In particular, there are at 
least two TeV blazars with redshift similar (PKS 2155-304 at z = 0.1 17) or higher (1ES 2344+514 
at z = 0.44) than H1426+428; but the data which are available for these blazars (Roberts et al. 
1998, Catanese et al. 1998) are not sufficiently accurate for our analysis. Finally, we need also 
new SPITZER data, to improve our knowledge of the Mid and Far Infrared background. 
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